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, ^ , Abstract. Valuable models for immortal solutions of Ricci flow that collapse 

■ with bounded curvature come from locally C7-invariant solutions on bundles 

Q N <—* Ai — ^->£5™, with Q a nilpotent Lie group. In this paper, we establish 
convergence and asymptotic stability, modulo smooth finite-dimensional center 
manifolds, of certain R^-invariant model solutions. In case N + n = 3, our 
results are relevant to work of Lott classifying the asymptotic behavior of all 
3-dimensional Ricci flow solutions whose sectional curvatures and diameters 
are respectively C?(f -1 ) and ©(t 1 / 2 ) as t — * oo. 
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1. Introduction 

There are many interesting open questions regarding the geometric and analytic 
properties of immortal solutions (M.,g(t)) of Ricci flow that collapse. One does 
not expect such solutions to converge smoothly in a naive sense. Instead, one may 
expect (A4,g(t)) to exhibit Gromov-Hausdorff convergence to a lower-dimensional 
geometric object. (See [llj . for example.) An impelling motivation for study- 
ing such solutions when the dimension of the total space is three is to obtain a 
better understanding of the role of smooth collapse in the Ricci flow approach to 
geometrization. See [19] and [14] for a broader discussion of some geometric and 
analytic problems related to collapse and their motivations. 

In a recent paper 20J, John Lott makes significant progress in understanding 
the long-time behavior of Ricci flow on a compact 3-manifold. He proves that any 
solution (M. 3 ,g(t)) of Ricci flow that exists for alH > with sectional curvatures 
that are and diameter that is ©(i 1 / 2 ) as t — > oo converges (after pull- 

back to the universal cover and modification by diffcomorphisms) to an expanding 
homogeneous soliton. A key component in the proof is the analysis of locally Q- 
invariant solutions of Ricci flow, where Q is a connected iV-dimensional Abelian Lie 
group. Such solutions are expected to constitute valuable models for the behavior 
of Type-Ill (immortal) solutions of Ricci flow in all dimensions, especially solutions 
undergoing collapse. We now review the basic setup from 20 . 

Let R N «-> M B n be a fiber bundle over a connected oriented compact base 
B, and let R N <-* £ -^-> B n be a flat vector bundle over B. Let £ x g M denote the 
fiber space \J beB {£b X M-b)- Assume that there exists a smooth map £ x e M. — > M 
such that over each basepoint b £ B, the map £b X A4 b — * M-b is a free transitive 
action. Assume also that this action is consistent with the flat connection on £ in 
the sense that for any open set U C B small enough that £\u ~ U X M. N is a local 
trivialization, the set 7r _1 (W) has a free M. N action, hence is the total space of a 
principal M. N bundle over U. In this way, M. may be regarded as an R^-principal 
bundle over B twisted by the flat vector bundle £. One can define a connection A 
on M. with the property that any restriction A\ v -\^ is an R^-valued connection. 

Let U C B be an open set small enough such that £\u ~ U x R N is a local 
trivialization and such that there exist a local parameterization p : W l — > U and a 
local section a :U — » 7r _1 (W). Then any choice of basis (ei, . . . , ejy) for R N yields 
local coordinates on 7r _1 (W) via 

M n x R N 3 (x a ,x l ) ^ (x l e t ) ■ a(p(x a )), 

where {x Q }2 = iare coordinates on 14, {x l }^ 1 are coordinates on W N , and • denotes 
the free K^-action described above. 

Let g denote a Riemannian metric on A4 with the property that this action is a 
local isometry. Then with respect to these coordinates, one may write 

n N 

(1.1) g= Yl 9 a pdx a dx +J2 G^idx* +Y / l =1 Al,dx a )(dx> +J2 n p =1 4dxP). 

q,/3— 1 id = 1 

Observe here that for b £ U C £>, ^™ « =1 g a p(b) dx a dx^ is the local expression of 
a Riemannian metric on B, A l a (b) dx a is locally the pullback by a of a connection 
on 7T gives a Euclidean inner product on the 

fiber over b. 
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A one-parameter family (Ai,g(t) : t £ T) of such metrics evolving by Ricci flow 
in a nonempty time interval I constitutes a locally M. N -invariant Ricci flow solution. 
In [20], Lott shows that such a solution is equivalent, modulo diffeomorphisms of 
M. and modifications of A by exact forms, to the systemQ 

(1.2a) J^ a/3 = -2R aP + ^& k & i V a G ij V (3 G ki + g^G^dA^dA)^, 
(1.2b) ^4 = -(6dAY a + &^G ]k {dA) k al 

(1.2c) ^Ga = AGy - G ke V a G lk V a G i3 - ^g ai 9 05 G lk G 3l {dA) k al3 {dA)i fS . 

Abusing notation, we denote a solution of (|1.2|) by g(i) = (g(i),A(t),G(t)). In 
order to study the long-time behavior of such systems, we engineer a transformation 
that turns certain model solutions into fixed points whose asymptotic stability can 
be investigated. To facilitate this, we assume that M. admits a flat connection, 
allowing us to regard A as an R^-valued 1-form. Let a function s(t) and a constant 
c be given. Let o~(t) be any positive antiderivative of s. Consider the transformation 

(g(t),A(t),G(t)) = g(t) » g(r) = (ff(r), A(r), G(r)), 

where 

g = a-~ x g, A = a-^A, G = a c G, t= a^^di, (t el). 

J t 

Observe that the exponent —(1 + c)/2 above is necessary so that no factors of a 
appear in the transformed system (|1.3|) below. Examples [TH3] below illustrate why 
this is desirable. Under the transformation g(£) i— > g(r), system (|1.2I) becomes^ 

(1.3a) 

= -2i? Q/3 + ^G ik & e V a G l3 VpG ke + g^G^idAY^idA)^ - sg aP , 

(1.3b) 

= -(*dA)j. + G^V^,- fe (dA)^ - ^sAi, 

(1.3c) 

O 1 

— G id = AGy - G M V a Gi k V a Gej - -g a ~' ' g ps G lk G 3l (dA)%(dA)\ s + csG^. 

We call this system a rescaled locally R. N -invariant Ricci flow. 

Example 1. On a solution (l w xB, g(i)) of (jl.2| that is a Riemannian product over 
a nonflat Einstein base (B, g), one may choose coordinates so that G is constant, 
A vanishes, and g(t) = — xtg(-x), where x = ±1 is the Einstein constant such 
that 2Rc(g(— x)) — xg(—x). (Note that g(t) exists for t < if x = 1 and for 
t > if x = — 1.) The choices s = — x, c = 0, and to = — x transform g(io) into a 
stationary solution g(0) of the autonomous system (|1.3p . 

iHere V a = V« = S^gf?, (dA)^ = V«-*j, - V^, (i<U}'„ = -V^dA)^, and 

AGy = g^VcVpGij = 9 af) { a jt dx ,i Gij - g^Gy), where f represents the Levi-Civita 
connection of g. 

2 Here V a = V« = 9^^, (dA)*^ = V Q A^ - VpA*,, (MA)* = - V 3 (dA)^ , and 

AGij = j^VaV^Gy = 9 ali ( a jt axl ] Gij - r^gfyGij), where T represents the Levi-Civita 
connection of g. 
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Example 2. A somewhat more general normalization is as follows. Suppose that 
g(r) may be regarded as a well-defined metric on B. Let V(r) denote the volume 
of (B,g(r)) and define 

r = R-\\VG\ 2 -\\dA\\ 

where everything is computed with respect to g. Because = — J B r dfi — ^sV(t), 
it follows that V is fixed if and only ifjfl 

s = (f r d/i . 

n Jb 

Consider a Riemannian product solution over an Einstein base. Such a solution 
may be written as g(t) = {git), A, G), where A vanishes identically and G is fixed 
in space and time. For any to in its time domain of existence, taking c = makes 
g(0) = (er -1 (£0)5(^0), 0, G) into a stationary solution of (| 1 . 3[) for any choice of 
a-^to) > 0. 

Example 3. The mapping torus M f of / : R N -> R N is ([0, 1] x R N )/ ~, where 
(0,p) ~ (1, f(p))- With respect to coordinates x° 6 S 1 and (a; 1 , . . . ,x ) € R N on 
the circle bundle R w ^ Mf — ^iS 1 , a locally R^-invariant metric has the form 

g = u dx° <g> da; + A t (dx° ® dec 1 + dec* ® da; ) + Gy da; 1 ® da; 3 , 

where u, A\ and G depend only on a; , and Ai = GijA^ . System (jl.3|) reduces to 

(1.4a) A« = \ \d G\ 2 - su, 

or 2 

d 1 + c 

(1.4b) ^ = -^sA*, 

(1.4c) ^-Gij = u-\vlG tl - G ke d G tk d G e3 ) + csG l3 , 

where |c?oG| 2 = G lk G^ l d G t jd Gk(. and V 2 , = <9q — r°, 9n. A motivating case 
occurs when N = 2. Compact 3-manifolds with sol geometry are mapping tori of 
automorphisms T 2 — > T 2 induced by A £ SL(2, Z) with eigenvalues A -1 < 1 < A. 
Now let g be a locally homogeneous sol-Gowdy metric of the kind considered as 
initial data in [IB] and |17j . so that A = 0. Without loss of generality, one may 
parameterize the initial data by arc length (thereby setting u = 1 at t = 0) and 
thereby write g = da; <g> dx° + Gij dx % da;- 7 , where one identifies G with the matrix 

r e F+w 

e F ~ w 

Here F is an arbitrary constant and 8qW = 2 log A is a topological constant deter- 
mined by the holonomy A around S 1 = [0, 1]/ ~. It is easy to check that this g 
becomes a stationary solution of (| 1 .4[) if and only if s — (2 log A) 2 > and c = 0. 
For the corresponding unsealed solution g = (u(t), A(t),G(t)), one has u = 1 + st, 
with A = and G fixed in time. 

Our purpose in this note is to investigate the stability of certain natural fixed 
points of system (ll.3[) like those in Examples H])-©. More precisely, we establish 
convergence and asymptotic stability of rescaled Ricci flow, modulo smooth finite- 
dimensional center manifolds, for all initial data sufficiently near certain model 



^Given any smooth function / : B — ► K, we define <fgfd[j,= ' •g~ ■ 
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solutions in the cases N = 1 or n = 1. These are the cases most amenable to study 
by linearization and also most relevant to applications when the total dimension is 
N + n = 3, e.g. to 20J . Because stability of flat metrics (modulo smooth finite- 
dimensional center manifolds) was established in [13| . we assume in what follows 
that the model solutions in question are not flat. 

A standard method for establishing stability of a nonlinear flow g r = F(g) near 
a stationary solution g is to proceed in two steps. (1) Compute the linearized 
operator Lg and establish that it is sectorial with stable spectrum Q (2) Deduce 
stability of the nonlinear flow in appropriate function spaces from properties of its 
linearization. For Ricci flow, this method is often not as easy as one might expect. 
There are several reasons for this. To apply existing theory such as [6], one requires 
the linearization to be elliptic. But because Ricci flow is invariant under the full 
infinite-dimensional diffeomorphism group, its linearization Lg is not elliptic unless 
one fixes a gauge, typically by the introduction of DeTurck diffeomorphisms. Even 
then, Lg may fail to be self adjoint. (For example, see [14] and Section [5] below.) 
Furthermore, in many cases of interest, such as those considered in 13J and in this 
paper, the spectrum of Lg has nonempty intersection with the imaginary axis. In 
general, this allows families of stationary or slowly changing solutions to compose 
local C r center manifolds, which can shrink as r increasesjj Finally, the DeTurck 
diffeomorphisms can introduce instabilities that may not be present in the original 
equation. (For example, the DeTurck diffeomorphisms solve a harmonic map flow. 
But the identity map id : S n — > S n of the round sphere is unstable under harmonic 
map flow for all n > 3. See Remark [2] below.) In spite of these obstacles, stability 
results are obtained in [13] and (by a somewhat different method) by Sesum in [23) . 

In this paper, we prove the following results that imply convergence in little- 
Holder spaces to be defined in Section I2TT1 below: 

Theorem 1. Let g = (<?, A, u) be a locally M 1 -invariant metric of the form \3.1\) on 
a product M 1 x B, where B is compact and orientable. Suppose that g has constant 
sectional curvature — l/2(n — 1), A vanishes, and u is constant. Then for any 
p G (0, 1), there exists 9 € (p, 1) such that the following holds. 

There exists a (1 + 9) little-Holder neighborhood U of g such that for all initial 
data g(0) G U, the unique solution g(r) of x-rescaled locally M 1 -invariant Ricci flow 
i3.3\) exists for all r > and converges exponentially fast in the (2 + p)- Holder norm 
to a limit metric goo — (500,^007^00) such that g x is hyperbolic, A^ vanishes, and 
Woo is constant. 

Theorem 2. Let g = (g,A,u) be a locally ^-invariant metric of the form f^. 1\ ) 
on a product M 1 x S 2 . Suppose that g has constant positive sectional curvature, A 
vanishes, and u is constant. Then for any p G (0,1), there exists 9 G (p, 1) such 
that the following holds. 

There exists a (1 + 9) little-Holder neighborhood U of g such that for all initial 
data g(0) G U , the unique solution g(r) of volume-rescaled locally ^-invariant 
Ricci flow H-'ity exists for all t > and converges exponentially fast in the (2 + p) - 
Holder norm to a liinit metric goo 

= (<7oo,Ax>, Woo) such that g^ has constant 
positive sectional curvature, ^4oo vanishes, and Uqo is constant. 

4 See Definition [T] on page HOl 

^For example, the "inner layer" asymptotics of a Ricci flow neckpinch correspond to non- 
stationary solutions in the kernel of its linearization at the cylinder soliton. See [l] . 
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Theorem 3. Let g = (u, A, G) be a metric of the form 15. 1\) on the mapping torus 
R N ^ Ma^S 1 of a given nonzero A G Gl(iV, R) admitting a flat connection. 
Suppose that u — 1, A vanishes, and G is a harmonic- Einstein metric i5.5\) . Then 
for any p G (0, 1), there exists 9 G (p, 1) such that the following holds. 

There exists a (1 + 6) little-Holder neighborhood U of g such that for all initial 
data g(0) G U , the unique solution g(r) of holonomy-rescaled locally M. N -invariant 
Ricci flow i f 5. £j) exists for all t > and converges exponentially fast in the (2 + p)- 
Holder TiOTTTi to a limit uictric goo 

= (uoo , Aoo , Goo ) such thatuoo = 1. A^ vanishes, 

and Goo is harmonic- Einstein. 

This paper is organized as follows. In Section [21 we review a general theory of 
asymptotic stability for quasilinear pde in the presence of center manifolds, then 
establish the context in which that theory applies here. In Sections [31 HI and0 we 
study and prove stability of >r-rescaled flows, volume-rcscalcd flows, and holonomy- 
rescaled flows, respectively. Each of these is a suitably-chosen variant of (|1.3jl . In 
Remarks [3] and HI below, we provide (partial) explanations of these apparently ad 
hoc choices of normalization. 

Acknowledgement. I warmly thank John Lott for sharing his results with me 
and raising the questions that are addressed in this note. 

Acknowledgement. I also thank James Isenberg for introducing me to several 
problems closely related to those considered in this note. Motivated by Kaluza- 
Klein theories in physics, Dedrickson and Isenberg did early work on locally K 1 - 
invariant solutions (both with and without a connection), as well as generalizations 
of Example [3] to higher-dimensional fibers (without a connection). 

Acknowledgement. My sincere thanks go to the referees for their suggestions to 
improve the exposition and notation of this paper. 

2. Center manifold stability theory for quasilinear systems 

In this section, we recall relevant aspects of the theory that allows one to derive 
rigorous conclusions about the asymptotic behavior of the nonlinear system (]1.3|) 
from its linearization. Good sources are [6] and [25]. Our approach here mainly 
follows the latter. The reader familiar with [13] can safely skim this section. 

2.1. Maximal regularity spaces. Because optimal asymptotic stability results 
are obtained in continuous interpolation spaces, we begin by recalling aspects of 
the maximal regularity theory of Da Prato and Grisvard [5] . (Section 12.21 below 
provides some motivation for our use of this theory.) 

Let y be a compact Riemannian manifold, possible with boundary dy. The main 
cases we have in mind in this note are: (1) 3^ is a compact hyperbolic manifold; (2) 
y is the round sphere; or (3) y = [0, 1]. 

Let T, Q (y), T, 1 (y), and E 2 (3^) respectively denote the spaces of C°° functions, 
C°° 1-forms, and C°° symmetric (2,0)-tensor fields supported on y. If dy ^ 0, 
as in case (3), we restrict to those functions, forms, and tensor fields satisfying 
prescribed linear boundary conditions like (|5.8|) - (|5.10p . respectively. 

Given r G N and p G (0,1), let Z r +P (y), ^"{y), and Y? 2 +P {y) denote the 
closures of T, (y), Si (y), and T, 2 (y), respectively, with respect to the relevant 
r + p Holder norms. These are the well-known little-Holder spaces. (Recall that 
C°° representatives are not dense in the usual Holder spaces.) Any little-Holder 
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space \f +p is a Banach space. Moreover, if s < r and a < p, then there is a 
continuous and dense inclusion fj r+p <—* fj s+cr . 

Given a continuous dense inclusion B\ Bo of Banach spaces and 9 G (0, 1), 
the continuous interpolation space {Bo, B\)g is the set of all igBj such that there 
exist sequences {y n } C Bq and {z n } C -Bi satisfying 



(2-1) ^ 



T" 

|j/n|ls = o(2-" e ) asn^oo 



||z„|| Bi =o(2"( 1 - 9 )) as^w. 

Continuous interpolation spaces are introduced in [5]. By [9], they are equivalent 
in norm to the real interpolation spaces also found in the literature. In fact, the 
norm \\x\\@ on (Bo,B\)g is equivalent to 

inf (sup{2" 9 \\y n \\ , 2 -"( 1 ~ e ) \\z n \\ } : {y n } and {z n } satisfy 

ln£N 

The little-Holder spaces are particularly well adapted to continuous interpola- 
tion. (See [27].) Indeed, let () r+p and t) s+a be little-Holder spaces. Let s < r G N, 
< cr < p < 1, and < 9 < 1. If 0(r + p) + (1 - 0)(s + <r) ^ N, then there is a 
Banach space isomorphism 

(2.2) (f) s+CT , ff +p )e = ()[ 0r +( 1 - 9 ) s ]+[ 9 P+( 1 - e ) ,T ] 

and there exists C < oo such that for all 77 G \) r+p , one has 

(2-3) IMI(i>»+- 1 s,'-+ P ) e ^ C ll r /llj~+^ h\\\,r+P ■ 

These properties make the little-Holder spaces highly useful for our purposes. 

2.2. A general center manifold theorem. To prove the results in this note, we 
invoke a special case of a theorem of Simonett 25j, whose hypotheses we now recall. 
(Compare 13J . ) Our notation is as follows. If X and Y are Banach spaces, then 
£(X, Y) is the set of bounded linear maps X — > Y. If (X, d) is a metric space, then 
-B(X, x, r) is the open ball of radius r > centered at x G X. If L is a linear operator 
on a real space, we denote its natural complexification by L c (m + iv) = Lu + iLv. 

If X ^ Y is a continuous and dense inclusion of Banach spaces and Q : X — > Y 
is a nonlinear differential operator, X9gn Q(g) G Y, we denote its linearization 
at g G X by Lg = Q'(g) : ©(Lg) CY^Y. Our main assumptions below are that 
g 1 ► Q(g) is a quasilinear differential operator and that Lg : g 1— > Lg(g) generates 
an analytic strongly-continuous semigroup on £(Y, Y). 

To apply Simonett 's stability theorem, one must verify certain technical hypothe- 
ses with respect to a given quasilinear parabolic equation and ordered Banach spaces 
satisfying 

(2.4) Xi C Ei G X G E and X x C X Q G X^ G X , 

whose precise relationships we describe below. The reasons for introducing this 
curious arrangement of Banach spaces come from the beautiful maximal regularity 
construction of Da Prato and Grisvard [5] . As motivation, consider a linear initial 
value problem 

r u(t) = L U (t) + f(t), 

[ ) \ u(0)=uo, 
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posed on a Banach space Yo, where the linear operator L generates a strongly- 
continuous analytic semigroup on £(Yo,Yo). Then L is a densely-defined closed 
operator whose domain, when equipped with the graph normQ naturally becomes 
a Banach space, which we may call Yi. Thus Yi > Yo is a continuous and dense 
inclusion. Now (for some fixed T > 0) suppose that / : [0, T] — > Yo is a bounded 
continuous function. Then the formal solution of (*) is given by the integral formula 



However, the convolution term above may not have enough regularity to make this 
formal solution rigorous unless one imposes stronger hypotheses (like requiring, 
for example, that / : [0, T] — > Yi be bounded and continuous). Da Prato and 
Grisvard's theory overcomes this difficulty by restriction to suitably chosen inter- 
polation spaces. In these (necessarily non-reflexive) spaces, ii and Lu enjoy the 
same regularity as /, which justifies the language "maximal regularity" . 
Stability theory for nonlinear autonomous initial value problems 



near a stationary solution N(0) — often proceeds by linearization, in which one 
replaces (**) by (*) with L = N'(0) and f(t) = N(u(t)) - Lu(t). Here the regu- 
larity of / is of evident importance. This is where one exploits the (|2.4|) hierarchy. 
Suppose that (the complexification of) L : Ei — > Eo is a sectorial operator^ The 
arguments Simonett uses to prove Theorem U below require the linear problem (*) 
to have a bounded continuous solution u : [0, T] — ► Xi for every bounded continuous 
/ : [0, T] — > Xo- This will fail for general sectorial Ei — > Eo but can be achieved by 
restricting L : Xi — > Xo to well-chosen continuous interpolation spaces Xi °-> Xo. 

A similar stability theorem is derived by Da Prato and Lunardi [6]. We use 
Simonett 's theorem because it takes optimal advantage of the parabolic smoothing 
properties of the quasilinear equation (|2.5[) in continuous interpolation spaces [25j . 
With regard to (|2.4[) . these properties allow one to show that invariant manifolds 
can be exponentially attractive in the norm of the smaller space Xi for solutions 
whose initial data are close to a fixed point in the larger interpolation space X a . In 
particular, such solutions immediately regularize and belong to Xi for all t > 0. 

With this brief introduction in hand, we are ready to state the seven hypotheses 
one needs to apply the theorem. In our applications, we choose certain little- Holder 
spaces whose properties, e.g. (|2.2|) . greatly simplify verification of the hypotheses. 
While perusing those hypotheses, readers may wish to consult (|2.6p and (|2.7|) below 
for specifics of how the various spaces are realized in the remainder of this paper. 

Assumption 1. Xi Xo and Ei Eo are continuous dense inclusions of Banach 
spaces. For fixed < (3 < a < 1, X a and "Kp are continuous interpolation spaces 
corresponding to the inclusion Xi °-> Xo- 

Assumption 2. There is an autonomous quasilinear parabolic equation 



'The graph norm with respect to Yq of a suitable linear operator L : D(L) C Yq — > ¥q is 




< t < T. 




(2.5) 



«-g(r)= Q(I(r)), 



(r > 0), 



X 



\d(l) — IfIIyq + II^IIyq' 

^See Definition [l] on page HOl 
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with the property that there exists a positive integer k such that for all g in some 
open set G/3 C X/3, the domain D(Lg) of Lg contains Xi and the map g i— > Lg|x a 
belongs to C k {G,3,£(K u X )). 

Assumption 3. For each g G &p, there exists an extension Lg o/Lg to a domain 
0(g) that contains Ei (hence is dense in Eq^. 

Assumption 4. For eac/i g £ G a = G^ nX , L§|ei G £(Ei,Eo) generates a 
strongly- continuous analytic semigroup on £(Eo,Eo). (Observe that for g G G Q , 
this implies that ©(g) becomes a Banach space when equipped with the graph norm 
with respect to Eo .) 

Assumption 5. For each g G G Q , Lg is the part o/Lg in XoB 

Assumption 6. For each g G G Q} there exists 8 G (0, 1) such that Xo = (Eo, B)(g))e 
and X! S (Eo,B(g))i+fl, toAere (E ,B(g))i+ e = {g G D(g) : Lg(g) G (E ,D(g)) e } 
as a set endowed with the graph norm o/Lg to£/i respect to (Eo,D(g))g@ 

Assumption 7. Ei ^ Eo is a continuous and dense inclusion such that 

there exist C > and (5 G (0, 1) suc/i t/iat /or a^/ ?y G Ei, one has 

Simonett obtains (a more general version of) the following result: 

Theorem 4 (Simonett). Let Lg denote the complexification of the linearization 
L g of \2.5}) at a stationary solution g of h2.5\) W\ Suppose there exists A s > such 
that the spectrum a of Lg admit the decomposition a — o~ s U {0} ; where is an 
eigenvalue of finite multiplicity and a s C {z : Rez < — A s }. If Assumptions H-0 
hold, then: 

(1) For each a G [0,1], there is a direct-sum decomposition X Q = X^ © X^j, 
where X^, is i/ie finite- dimensional algebraic eigenspace corresponding to 
the null eigenvalue of Lg . 

(2) For each r G N, there exists d r > suc/i £/iai /or aZZ d G (0, d r ], i/iere 
exists a bounded C r map 7^ : _B(X^,g, d) — > X^ such that 7 J(g) = and 
Z?72(g) = 0. The image of 7^ /ies in i/ie closed ball _B(X^,g, d). Its graph 
is a /oca? C r center manifold r[ oc = {(/i, : ^ G -B(^i,g, d)} C Xi 
satisfying T g T^ oc = XJ. Moreover, r}" oc is invariant for solutions of i2.5\) 
as long as they remain in _B(X^,g, d) x _B(X^,0, d). 

(3) Fke A G (0, A s ). Then for each a G (0, 1), there exist C > and d G (0, d r ] 
such that for each initial datum g(0) G S(X Q , g, d) and aZZ times r > smc/i 
£/ia£ g(r) G B(K a , g, d), i/ie center manifold T\ oc is exponentially attractive 
in the stronger space Xi in the sense that 

|k s g(r) - 75(^ c g(r))|| Xl < -fee" AT lk s g(0) - 75(^1(0)) || Xq • 



If X is a Banach space with subspace Y and L : D(L) C X — > X is linear, then L Y , t/ie part of 
L in Y, is defined by the action L Y : x 1— » on the domain D(L Y ) = {x £ D(L) : Lx £ Y}. 

^To see that this assumption is equivalent to condition (n) in 1251 Section 4] , recall the fact (see 
e.g. [2 1 1 Propositions 2.2.2 and 2.2.4]) that if A is the generator of a Co analytic semigroup in a 
Banach space X and < 9 < 1, then (X, HJ(A)) 9 S »a(9) = {x £ X : lim txo t~ e (e tA x - x) = 0}. 

^Note that L g is the operator that appears in Assumption 2. 
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Here, g(r) is the unique solution of i2.5\) . while ir s and 7r c denote the pro- 
jections onto = (X^,Xq) q and X°, respectively. 

2.3. Prerequisites for application of the theorem. We now establish the con- 
text in which we apply Theorem 21 

Let g T = Q(g) denote the rescaled IR^-invariant Ricci flow system (|1.3|) . modified 
by DeTurck diffeomorphisms with respect to a background metric g chosen as in 
Sections [3H5]below. We assume that g is a smooth stationary solution, namely that 

gr = Q(g) = 0. 

For fixed < a < p < 1, consider the following nested spaces: 

T (y) 

u 

x = x° +p (y) x s^oo x 

(2.6) 



E 


— ^0 


- a {y) x y a + a {y) x 


U 






X 


- ^0 


vp (y) x TP+ p {y) x e°+'i 


u 






Ei 


— ^0 


- a {y) x y 2+a {y) x y 2 + a 


U 






Xi 


— ^0 


hp (y) x y\ +p {y) x e^+'i 



For fixed l/2</3<a<l and e > to be chosen, define 

(2.7) X^ = (X ,Xi)^ and X a = p ,X l ) a 
and 

(2.8) Gp = B(X fj ,g,e) and G„=G^nX a . 

For each g G G/3, let Lg denote the linearization of Q at g, regarded as an 
unbounded linear operator on Eo with dense domain D(g) = Ej. Let Lg denote 
the corresponding operator on Xo with dense domain B(Lg) = Xi. 

Recall the following: 

Definition 1. A densely-defined linear operator L on a complex Banach space 
X ^ {0} is said to be sectorial if there exist a £ (tt, 2tt), lu G R, and C > such 
that the "sector" 

(2.9) S Q , W = {A e C : A ^ w, |arg(A - w)| < a} 
is contained in the resolvent set p(L), and such that 

holds for all A £ S a ,u Q p(L). 

The following lemmas verify the technical hypotheses needed for Theorem |4j 
Lemma 1. 

(1) g 1 — > Lg is an analytic map Gfj — > £(Xi,Xo). 

(2) g 1— > Lg is an analytic map G Q — > £(Ei,Eo). 

(3) // L^ is sectorial, then there exists e > suc/i £/ia£ /or a/Z g m f/ie set 
G Q defined by i2.8\) , L? is ifte infinitesimal generator of an analytic Co- 
semigroup on £(Eq,Eo). 
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Proof. Statements (l)-(2) are proved in Lemma 3.3 of [13]. 

To prove statement (3), first observe that Lg generates an analytic Co-semigroup 
on £(Eq,Eo): it is a standard fact that a sectorial operator generates an analytic 
semigroup; one knows that semigroup is strongly continuous because Lg is densely 
defined by construction. Now if g € G a , then by statement (2), we can choose 
e > small enough so that 

|£g-££ 

where C > is the constant in (|2.10p corresponding to Lg. As is well known, 
this implies that L? is sectorial as well. (See [21] Proposition 2.4.2], for example.) 
Statement (3) follows easily. □ 

Lemma 2. If Lg is sectorial, then the choices made in S2.6]) - [278\) ensure that 
Assumptions^^ are satisfied for the system g T = Q(g) given in (1.3\) . 

Proof. Assumption 1 holds by construction. Assumption 2 follows directly from 
statement (2) of Lemma [TJ because the "off-diagonal terms" in (| 1 .3|) are all con- 
tractions involving at most one derivative of g. Assumptions 3 and 5 hold by 
construction of L g and L g . Assumption 0] is a consequence of the hypothesis on Lg 
and statement (3) of Lemma [T] To verify Assumption 6, first observe that putting 
6 = (p — a) /2 G (0, 1) into the isomorphism identity (|2.2p yields Xo = (Eo,Ei) e ; 
then recall that standard Schauder theory implies that the graph norm of Lg with 
respect to Eo is equivalent to the Ei norm. (So the graph norm of L g with respect 
to Xq is equivalent to the Xi norm.) Finally, Assumption 7 is implied by (|2.3p . □ 



< 77- 



3. The case N = 1: ^-rescaled locally ^-invariant Ricci flow 

To consider evolving Riemannian metrics on a bundle R 1 A4-^>B n , it is 
convenient to change notation. Let (x , ... , x n ) be local coordinates on B n and 
let x° denote the coordinate on an R 1 fiber. If (jM,g(t) : t G T) is a locally R 1 - 
invariant solution then there exist a Riemannian metric g, a 1-form A, and a 

function u, all defined on B n x T, such that we can write g in local coordinates as 

(3.1) g = e 2a dx ®dx + e 2a A k (dx° '®dx k + dx k ®dx ) + (e 2a ■A i A j +g ij ) dx i ®dx° . 

We again abuse notation by writing g(f) = (g(t), A(t),u(t)). In the coordinates of 
this section, system (|1.2[) takes the form 

d 

(3.2a) — g i:j = -2R tJ + e 2u g M (dA) lk (dA) je + 2V.it Vj-u, 

(3.2b) ^A k = g^{Vi(dA) jk + 2(Viu)(dA) jk }, 

(3.2c) = Au-^e 2a \dA\ 2 . 

We are interested in the stability of metrics of this type on trivial (i.e. product) 
bundles. It turns out that the most practical rescaling depends on the sign of the 
curvatures of (B,g). We explain this in Remarks [3HH below. 
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3.1. >r-rescaled flow. Suppose that (_M,g) is a Riemannian product over an Ein- 
stein manifold with 2Rc(g(— x)) = Kg{—x) for some x = ±1. Then the renormal- 
ization in Example [1] yields 



(3.3a) 
(3.3b) 
(3.3c) 



d 

—g l3 = -2R i3 + e 2u g M {dA) lk {dA) ]t + 2V,uV 3 u + x 9lj , 
— A k = g»{Vi{dA) jk + 2(V<u)(dA) ifc } + ^A k , 



9 

— — ii = Au 

OT 



\e 2u \dA\\ 



We call this system (in which all quantities are computed with respect to g) the 
x-rescaled locally R 1 -invariant Ricci flow. It is most satisfactory for studying the 
case x = — 1. (In Section 01 we apply a different rescaling that is more suitable for 
the case h = +1.) 



3.2. Linearization at a stationary solution of >f-rescaled flow. Any Rie- 
mannian product (R x B,e 2u dx° <g> dx° + g^ dx % <8> da;- 7 ') with 2Rc(g) = xg, A 
identically zero, and u constant in space is clearly a stationary solution of (|3.3[) . 
Let {g + h, B,u + v) denote a perturbation of such a fixed point g = (<?, 0, u). De- 
fine H — tr g h, and let Ae denote the Lichnerowicz Laplacian acting on symmetric 
(2, 0)-tensor fields. In coordinates, 



(3.4) 



Afhj 



Ah, 



2-^ipqjh 



Lemma 3. The linearization of \3. 3\) at a fixed point g = (g, 0, u) with 2 Rc = Kg 
and u constant acts on {h, B, v) by 



(3.5a) 
(3.5b) 
(3.5c) 



dr ij 



Aghij + {Vi^/i)^ + V-j{8h)i + ViVjH} + xh tj , 



—B k = -(SdB) k + ^B k , 
d 

—v = Av. 

OT 



Proof. For the convenience of the reader, we begin by recalling a few classical 
variation formulas (all of which may be found in [21 Chapter 1, Section K], for 
instance). Let g(e) be a smooth one-parameter family of Riemannian metrics such 
that 



9(P) = 9 



A 9 

and — 

Oe 



fJ 



h. 



e=0 



Using tildes to denote geometric quantities associated to g and undecorated char- 
acters to denote quantities associated to g, one computes that 



0_ 

ds 
d_ 

Ik 



f)' J 



-h 13 



-g lk g 3l hu) 



e=0 



(V^ + V,/if-V%), 



e=0 



i? 



£ = 



\7i(8h)j + Vj{Sh)i + V.VjH} . 



The lemma follows by careful but straightforward application of these formulas. □ 
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The linear system p.5[) is autonomous but not quite parabolic. We impose 
parabolicity by the DeTurck trick [8] . Fix a background connection T and define 
a 1-parameter family of vector fields W(r) along a solution g(r) of (|3 . 3|) by 

(3.6a) W° = SA, 

(3.6b) w k -/'(it r;; ;i (l<fc<n). 

The solution of the x-rescaled K -invariant Ricci-DeTurck flow is the 1-parameter 
family of metrics 7/>*g(r), where the diffeomorphisms ip T are generated by W(t), 
subject to the initial condition -0o = id. In what follows, we take T to be the Levi- 
Civita connection of the stationary solution around which we linearize. Observe 
that a stationary solution g = (g,0,u) of (|3.3j) with 2Rc = xg and u constant is 
then also a stationary solution of the ^-rescaled Ricci-DeTurck flow. 

Lemma 4. The linearization of the x-rescaled Ricci-DeTurck flow at a fixed point 
g = (<?, 0,u) with 2Rc = xg is the autonomous, self-adjoint, strictly parabolic 
system 

d 

— (h B v) = L (h B v) = (L 2 h LiB L v) , 

where 

(3.7a) L2/1 = Aih + xh, 

(3.7b) L 1 B = A 1 B+^B, 

(3.7c) L v = A v. 

Here — A p = dS + 5d denotes the Hodge-de Rham Laplacian acting on p-forms. 

Proof. The Christoffcl symbols T for a fixed point g = (g, 0, u) vanish if any index 
is zero and satisfy = r*j- otherwise, where 1^ are the Christoffel symbols for g. 
Hence the DeTurck correction terms to the linearization (|3.5I) of (|3 . 3|) are 

(£ w g)oo = 0, 
(Avg)ofc = e 2u (d6B) k , 
{Cwsh = Vi(W)j + VjiSfyi + ViVjH, 
whence the result is immediate from (13.11) and (13.51). □ 



3.3. Linear stability of >f-rescaled flow. Now let us make the stronger assump- 
tion that g = (3, 0, u) is a fixed point of the >r-rescaled R 1 -invariant Ricci-DeTurck 
flow with u constant and g a metric of constant sectional curvature — 1/2 (n— 1). 

We say a linear operator L is weakly (strictly) stable if its spectrum is confined 
to the half plane Re z < (and is uniformly bounded away from the imaginary 
axis). 

Because L is diagonal, we can determine its stability by examining its component 
operators. The conclusions we obtain here will hold below when we extend L to a 
complex- valued operator on a larger domain in which smooth representatives are 
dense. 

Lemma 5. Let g = (g, 0, u) be a metric of the form \3.1\) such that u is constant 
and g has constant sectional curvature — l/2(n — 1). Then the linear system {3. 7| ) 
has the following stability properties: 

The operator Lq is weakly stable; constant functions form its null eigenspace. 
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The operator Li is strictly stable. 
If n > 3, then L 2 is strictly stable. 

If n = 2, then L 2 is weakly stable. On an orientable surface B of genus 7 > 2, the 
null eigenspace of Ai — 1 is the (67 — 6) -dimensional space of holomorphic quadratic 
differentials. 

Proof. The statements about Lo and Li are clear. 

Let Q = Q(h) denote the contraction Q = Rijkih ,h? and define a tensor field 
T = T(h) by 

(3.8) Tij k = Wkhij - Vih jk . 

Then integrating by parts and applying Koiso's Bochner formula [18] . one obtains 



(L 2 / l » = -||V^|| 2 + 2 / Q(h)d(x 

JB 



'" "Tf-\\Sh\\ 2 -l\\hf + / Q(h)dn 



2 



(3-9) = -~ \\T\\ 2 - \\Shf - \\Hf -JLJL \\h\\ 2 . 

2 2{n — 1) 2{n — 1) 

This proves that (|3.7a.|) is strictly stable for all n > 2. 

When n = 2, it follows from (|3.9p that a tensor field h belongs to the nullspace 
of L 2 = + R if and only if h is trace free and divergence free, with T(h) = 
vanishing identically. We first verify that the latter condition is superfluous. In 
normal coordinates at an arbitrary point p S B, all components of T(h) vanish 
except possibly Tn 2 = -T 2 \\ = d 2 hn - d\hi 2 and T 22 i = -T i22 = d\hii — d 2 h 2 i. 
But one has d\hu + d 2 h 2 \ = and d\h\ 2 + d 2 h 22 = because h is divergence free 
and d\hu + d\h 22 — and d 2 hu + d 2 h 22 = because h is trace free. It follows 
easily that T vanishes at p, hence everywhere. 

Now let He a trace-free and divergence-free tensor field on an orientable surface 
B of genus 7 > 2. To identify the nullspace of L 2 , it is most convenient to work in 
local complex coordinates. The trace-free condition implies that 

h = f(z) dz dz + f(z) dzdz 

for some function / : B — * C, while the divergence-free condition implies that / is 
holomorphic. Hence h is a holomorphic quadratic differential. If 1Z is a Riemann 
surface obtained from a closed surface by deleting p disjoint closed discs and q 
isolated points not on any disc, the Riemann-Roch theorem implies that the space 
of holomorphic quadratic differentials on 1Z has real dimension 6(7 — 1) + 3p + 2q. 
The result follows. □ 

When n = 2, the nullspace for the system is classically identified with the cotan- 
gent space to Teichmiiller space. 

3.4. Convergence and stability of x-rescaled flow. We now prove Theorem[T] 
To obtain asymptotic stability of xr-rescaled K 1 -invariant Ricci flow (|3.2[) with 
h = — 1 from its Ricci-DeTurck linearization (|3.7p , we shall freely use the theory 
reviewed in Section [2] and in particular the maximal regularity and interpolation 
spaces chosen in Section [ 



Proof of Theorem]^ Our notation is as follows, g = (g,A,u) is a locally K 1 - 
invariant metric (|3.ip on a product IR 1 x B, with B compact and orientable, such 
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that g has constant sectional curvature — l/2(n — 1), A vanishes, and u is constant. 
g(r) is the unique solution of ^-rescaled locally R 1 -invariant Ricci-DeTurck flow 
corresponding to an initial datum g(0). Then there are diffeomorphisms (p T with 
ipo = id such that the unique solution g(r) of ^-rescaled locally R 1 -invariant Ricci 
flow (|3. 3[) with initial datum g(0) is given by 



The proof consists of four steps. 

Step 1. We prove that (the complexification of) the linearization L = L g of 
xr-rescaled locally RMnvariant Ricci-DeTurck flow is sectorial. Observe that L is 
strictly elliptic and self adjoint. By Lemma El one may fix u > such that its 
spectrum cr(L c ) satisfies cr(L c )\{0} C (— oo, — uj). Standard Schauder theory then 
implies that L c is sectorial. (See Lemma 3.4 in [13] for a more detailed argument.) 

Step 2. By Step 1 and Lemma [21 we may apply Theorem 2] to g(r). Recall 
that X Q = (Xo,Xi) a , where Xo and Xi are defined by (|2.6p . Theorem [4] implies the 
following three statements: 

(1) For each a £ [0, 1], there is a direct sum decomposition X Q = W a © X^,, 
where X^j is the nullspace of L of dimension 



(Recall that if n = 2, then B must be an orientable surface of genus 7 > 2.) 

(2) For each r £ N, there exist d and a local C r center manifold r[ oc = 
graph (7S : B(Xf,g, d) — ► Xf ) which is invariant for solutions of Ricci- 
DeTurck flow as long as they remain in B(Xf , g, d) x B(K\ ,0, d). 

(3) For each a £ (0,1), there exist positive constants M — M(a) and d — 
d(a,r) such that for each initial datum g(0) £ £?(X Q ,g, d) and all times 
r > such that g(r) £ B(K. a , g, d), the center manifold r[ oc is exponentially 
attractive in the stronger space Xi in the sense that 



(3.10) \\n s g(r) - 7 5(T c g(r))|| Xl < —e~" T \\n s g(0) - YA^s(0))\ 



where 7r s and ir c are projections onto X s a = (X| , Xq) q and XJj, respectively. 

(All constants introduced here and below implicitly depend on g.) 

Step 3. We prove that the local center manifolds r[ oc coincide for all r, namely, 
that there is a unique smooth center manifold T = graph (7 : S(Xf , g, do) —> X|) 
consisting of fixed points of ^-rescaled locally R 1 -invariant Ricci flow (|3.3j) . First 
observe that for all n > 2 and all c £ R, (g, 0, c) is a stationary solution both 
of (|3.3p and Ricci-DeTurck flow. By (|3.10p . any such metric sufficiently near 
g = (g,0, u) must belong to all Fj r oc . When n = 2, Teichmiiller theory shows that 
there is a 6(7 — l)-dimensional space V of metrics g' = (g' , 0, u) near g such that 
g' is hyperbolic. Each such metric is a fixed point of ()3.3|) . hence evolves under 
Ricci-DeTurck flow only by diffeomorphisms. By diffeomorphism invariance, g'(r) 
remains hyperbolic, so g'(r) remains in V . By (|3.10p . g'(r) must belong to all rj" oc . 

Step 4. Fix a £ (1/2,1 — p/2). Then by the interpolation isomorphism (|2.2p . 
||-|L is equivalent to a (l + 0)-H61der norm, with 9 £ (p, 1). For all g(0) sufficiently 
near g, we now prove that g(r) converges in the X Q norm to an element of F. (We 
give the proof in general, even though certain steps are much easier when n > 2.) 



g(r) = #(g(r)). 




71/ 



r 
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Fix A 6 (0, u). Then by (f3TT0|> and Step 3 there exists C = C(M, A) such that 

(3.11) |k s g(r) - 7 (7r c g(r))|| Xi < Ce^ ||g(0) - g|j^ 

for all r > such that g(r) € B(K a ,g,d). Let S,e be positive constants to be 
determined so that < e < S < d, and suppose that g(0) € B(X Q ,g,e). Then it 
follows from (f3~TT]) and (O that 



(3.12) 



3 ~< s 



< C n Ce- AT e. 

Q+P 

If S/e is sufficiently large, this implies that 

l|g(r) ~ g|l Xl < e (l + ^f) < S, 

uniformly in time. Because g(r) = tp*(g(T)), the only way that the solution g(r) 
of Ricci-DeTurck flow could leave i?(X Q ,g, S) is by diffcomorphisms. But as was 
observed by Hamilton [15j . the diffcomorphisms tp T satisfy a harmonic map flow 

d 

—ip T = Ag (r)ig ^ r . 

with domain metric g(r) and codomain metric g. Because g(r) 6 _B(Xi,g,<5) for 
all r > 0, mild generalizations of standard estimates for harmonic map heat flow 
into negatively curved targets imply that the diffeomorphisms ip T exist for all r > 
and satisfy 

||Vr-id|| Xl < c 

for c = c(S). (See [TU] and [53], for instance.) Note that c does not increase if one 
makes 5 smaller. It follows that for 5 € (0, d) and e = e(5) £ (0, 8) sufficiently 
small, one has g(r) 6 B(X a , g, d) for all r > 0. By (|3.12j) . the result follows. □ 

4. The case N = 1: volume- rescaled locally ^-invariant Ricci flow 

In case (A4,g) is a Riemannian product over an Einstein manifold of positive 
Ricci curvature such that 

(4.1) g = e 2a dx°(E)dx + e 2i 'A k (dx°(g>dx k + dx k (g>dx ) + (e 2il A i A :j +g i:i )dx i (E)dx 3 , 

the rescaled ^-invariant Ricci flow system (|1.3[) corresponding to the normalization 
in Example [2] is a more suitable choice. 

4.1. Volume-rescaled flow. Volume-rescaled locally R 1 -invariant Ricci flow is 
the system 

2 f 

(4.2a) — 3ij = -2Rij + e 2u g M {dA) lk {dA) 3l + 2V,uV J u + -(6 r d/x ) gij , 
or n J B 

d . tSt „ „,„ 1 



(4.2b) — A k = gV{V l (dA) jk + 2(V l u)(dA) jk } + -(<b rd^)A k , 
or n J B 

(4.2c) A u = Au - -e 2u \dA\ 2 , 

or 4 

where all geometric quantities are computed with respect to g, and 
(4.3) r = R - ^e 2u \dA\ 2 - \Vu\ 2 . 
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4.2. Linearization at a stationary solution of volume-rescaled flow. Fix a 

background connection T and define a 1-parameter family of vector fields W(t) 
along a solution g(r) of (|4.2p by 

(4.4) W° = 5 A, 

W k =^{T%-t) (l<fe<n). 

The solution of the volume-rescaled M 1 -invariant Ricci-DeTurck flow is the 1- 
parameter family of metrics ^g(r), where the diffeomorphisms ip T are generated 
by W(t), subject to the initial condition ipo — id. Below, we take T to be the 
Levi-Civita connection of the metric about which we linearize. 

Any Riemannian product (Rx B, e 2u dx° ®dx° +gij dx l ®dx : > ) with g an Einstein 
metric, A identically zero, and u constant in space is clearly a stationary solution 
of the volume-rescaled M 1 -invariant Ricci-DeTurck flow. Let g = (g, 0, u) be such a 
fixed point, with Rc = Kg and hence r = nK. Again, let (g + h, B,u + v) denote a 
perturbation of g. We continue to write H — tr g h and to denote the Lichnerowicz 
Laplacian (|3.4[) by Af. Proceeding as we did for the xr-rescaled flow, one obtains 
the following. 

Lemma 6. The linearization of the volume-rescaled R 1 -invariant Ricci-DeTurck 
flow at a fixed point g = (g, 0, u) with Rc = Kg and u constant is the autonomous, 
self-adjoint, strictly parabolic system 
3 

— (h B v) = L (h B v) = (L 2 h LiB L v) , 

where 

(4.5a) L 2 h = A e h + 2K{h - -Hg\, 

n 

(4.5b) l, 1 B = A 1 B + KB, 

(4.5c) L v = A v. 

Here H — § B Hdu and — A p = dS + 5d denotes the Hodge-de Rham Laplacian 
acting on p-forms. 

Proof. The argument here needs a few more variation formulas in addition to those 
we recalled in the proof of Lemma [3] Once again, let g(e) be a smooth one- 
parameter family of Riemannian metrics such that 

d 

g(0) = g and — g = h. 

Using tildes to denote geometric quantities associated to g and undecorated char- 
acters to denote quantities associated to g, one recalls from |T3l Lemma 2.2] that 

R= -(AH - 5 2 h+ (Rc,h)), 

=o 



ds 
d_ 

Ik 

d 



dp, = -Hdfi, 



ds 



Rdfi = j>y^ \ R- j> Rdnj H - (Rc,/i)j dfj,. 

These formulas simplify nicely when Rc = Kg. Using this fact, the result follows 
by careful but straightforward calculation. □ 
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4.3. Linear stability of volume-rescaled flow. Now let us make the stronger 
assumption that g = (<?, 0, u) is a fixed point of the volume-rescaled K -invariant 
Ricci-DeTurck flow with u constant and g a metric of constant sectional curvature 
k > and hence K = (n — l)k. 

By passing to a covering space if necessary, we may assume that B n is the round 
n-sphere of radius y/l/k. Recall that the spectrum of the Laplacian Ao acting on 
scalar functions on (£>",<?) is {^j,k}j>o, where 

= ~jk(n + j - 1). 

The eigenspace Aj.k consists of the restriction to B n C IR n+1 of all j-homogeneous 
polynomials in the coordinate functions (x 1 , . . . , x n+1 ) of R n+1 . 

Given a symmetric (2, 0)-tensor field h, let / = f(h) denote its trace-free part 
defined by 

(4.6) h = -Hg + f. 

n 

It is then convenient to rewrite (|4.5aj) and (j4.5b|) as 

Tl — 1 — 

(4.7a) L 2 /i = Ah - 2kf + 2 k(H - H)g, 

n 

(4.7b) L 1 B = A 1 B+(n-l)kB, 

respectively!"! 

As above, the fact that L is diagonal lets us determine its stability by examining 
its component operators. The conclusions we obtain will continue to hold when we 
extend L to a larger domain in which smooth representatives are dense. 

Lemma 7. Let g = (,g,0, u) be a metric of the form 13.1]) such that u is constant 
and g has constant sectional curvature k > 0. Then the linear system |^.5[ ) has the 
following stability properties: 

Lo is weakly stable; its null eigenspace consists of constant functions. 

Li is strictly stable. 

If n = 2 , then L2 is weakly stable; its null eigenspace is 

{ipg : ip G A 0:fe u A ljfc ). 

If n > 3, then L2 is unstable. The sole unstable eigenvalue is (n — 2)k with 
eigenspace {tpg : ip G Ai^.}. The null eigenspace is {ipg : ip G Ao.fc} = {eg : c G R}. 

Proof. Weak stability of Lo is clear. 

It is well known (see [T^], for example) that the spectrum of Ax is 

{ A i,fc}j>i u i x 3,k - (n - 2)fc} J -> 1 . 

Because n > 2, it follows that the largest eigenvalue of Li is —k. 

Using (14. 6p . one obtains the decomposition L2/1 = -^LqH + L' 2 f where 

L' Q H = A H + 2(n - l)k(H - H), 

L' 2 f = A/ - 2kf. 

The operator L' 2 is evidently strictly stable. Because 

(L' a H,H) = - ||VF|| 2 + 2(n - l)k(\\H\\ 2 - VH 2 ), 



Recall that A is the rough Laplacian. 
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where V = Vol(B, g), one sees that the spectrum of L' is bounded from above by 
that of Ao + 2(n — l)k. Because A2,fc = — 2(n + l)k, the only possible unstable 
eigenfunctions for L' Q are elements of Ao^ and Ai < i~. It is easy to see that elements 
of Ao,fc (i.e. constant functions) belong to the nullspace of L' Q . And if H € Ai j, then 
one has H = and hence (L' H, H) = {Ai,* + 2(n - l)fc} \\H\\ 2 = (n - 2)fc ||H|| 2 . 
Because \\h\\ 2 = A ||ff|| 2 + ||/|| 2 , the result follows. □ 

Remark 1. The elements of {^3 : 95 6 Ai^} are infinitesimal conformal diffeomor- 
phisms (Mobius transformations): it is a standard fact that any h € {</?3 : </? G Ai.fc} 
satisfies 

Remark 2. The apparent instability of L2 for n > 3 is an accident of the DeTurck 
trick rather than an essential feature of Ricci flow. Indeed, Hamilton observed 
[15j that the DeTurck diffeomorphisms solve a harmonic map heat flow, with an 
evolving domain metric and a fixed target metric. It is well known that the identity 
map of round spheres S n is an unstable harmonic map for all n > 3. On the other 
hand, one can show by different methods that spherical space forms are attractive 
fixed points for normalized Ricci flow in any dimension. See [3] and [4j. 

In any case, the instability of L2 is in a sense a red herring, i.e. geometrically 
insignificant. Indeed, Moser |22j has shown that if g and g are any metrics on 
a compact manifold M such that Vol(JV, g) — Vol (A/ - , 3), then g = tp*g for some 
diffeomorphism i\) : M — » M and some metric g on M with djl = dfi. So, up 
to diffeomorphism, it suffices to consider perturbations that preserve the volume 
element d/i of a stationary solution g. These are exactly those h with H = 
pointwise. Restricted to such perturbations, L 2 is strictly stable. 

We can now at least partially explain our ad hoc choices of normalization. 

Remark 3. If g = (3,0, it) is a metric of the form (|3.ip such that u is constant 
and g has constant sectional curvature l/2(n — 1) > 0, then the operator L2 in the 
linearization (|3.7ll of >r-rescaled Ricci-DeTurck flow (|3.3j) satisfies 

(L 2 h, h) = i(|| Vff|| 2 + || 2 ) - IIV/II 2 - -J- ll/H 2 . 

So 1 is an eigenvalue with eigenspace {eg : c £ R}. 

Remark 4. If g = (3, 0, u) is a metric of the form (|3.ip such that u is constant and 
g has constant sectional curvature k < 0, then Koiso's Bochner formula [H] in the 
form 

\\Vh\\ 2 = l\\T\\ 2 + \\8h\\ 2 -nk\\ff 

implies that the operator L2 in the linearization (14. 5|) of the volume-rescaled Ricci- 
DeTurck flow gjg) satisfies 

(L 2 h, h) = -\ \\T\\ 2 - \\5h\\ 2 + (n- 2)k \\.f\\ 2 + 2 — k{\\H\\ 2 - VH 2 }. 
2 n 

Holder's inequality implies that the quantity in braces is nonnegative and vanishes 
exactly when \H\ is constant a.e. So L2 has a null eigenvalue with eigenspace 
{eg : c G K} in all dimensions. (The nullspace is of course larger when n = 2). 
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4.4. Convergence and stability of volume-rescaled flow. Here we exploit the 
fact that the identity map of the round 2-sphere is a weakly stable harmonic map. 

Proof of Theorem [H The proof is entirely analogous to the proof of Theorem [TJ 
except for Step 4. Here we must uniformly bound diffcomorphisms {ip T } solving 

d 

—<p T = A g(r)ig ^ T , 

^o = id:(5 2 ,g(0))^(5 2 ,g), 

where g is a metric of constant positive curvature and g(r) E -B(X l5 g,<5) for all 
t > 0. In this case, the required bound \\ip T — id|| Xi < c(8) follows from a result of 
Topping |26] . The remainder of the proof goes through without modification. □ 

5. The case n = 1: holonomy-rescaled locally R^-invariant Ricci 

FLOW 

We now consider evolving metrics on the mapping torus °-> M^-^S 1 of 
a given A 6 Gl(iV, R). We let x° denote the coordinate on S 1 « R/Z and take 
(ar, . . . , x ) as local coordinates on the fibers. Let c and s be constants and let V 
be a 1-form to be determined. Then the evolution of 

(5.1) g(r) = u da; <8> da; + A l (dx° ® da;* + da;* (g> da; ) + Gy dx l <8> dx J ' 

by rescaled locally R^-invariant Ricci flow (|1.3|) modified by diffeomorphisms gen- 
erated by the vector field V* is equivalent to the system 

d 1 2 
(5.2a) — m = - |<9 G| - su + (£ v »g) o, 

C»T 2 



(5.2b) = -^sA* + G^(C viS ) 3Q , 



(5.2c) ^-G VJ = u-^V^Gy - (9 G lfc )G fe£ (9 G J ,)} + caGy + (ZV»g)< 



We define 



(5.3) V = — (Sb«) <*c° + (Gy-a ^) da;* 



G, 
2 

where G > is a constant to be chosen below. Notice that V implements a DeTurck 
trick, as was done in Sections [3HU 

Hereafter, we drop the boldface font and write g for g. For later use, we note 
that the Christoffel symbols of g are 

1 

T 



(5.4a) 


1 00 


(5.4b) 


-pk 
1 00 


(5.4c) 


1 iO 


(5.4d) 


r°. 

y 


(5.4e) 


1 0j 


(5.4f) 





- ff 00 9 w + 5 0l 5o(G 4J ^), 



i ff feo a oU + .g fci ao(G l3 ^), 
^ G 4J , 



1 
2 

1 



^9 M d Q G 
1 
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5.1. Holonomy-rescaled flow. Following Lott [19] , we say G is a harmonic- 
Einstein metric if 

(5.5) Vldj = (d G tk )G ke (d G ej ). 

We henceforth assume that A4a admits a metric g of the form (|5.ip with G 
harmonic-Einstein. By an initial reparameterization, we may assume without loss 
of generality that u = 1 at t = 0. As we shall observe in Lemma [8] below, the 
harmonic-Einstein condition implies that |(?oG| 2 is constant in space. Motivated by 
Example[3l we choose constants c = and s = \ \QqG\ . Observe that the constant 
s is determined by the holonomy of A4\, because for any fiber metric G, one has 

(5.6) Ud = i ^A^Ap M \ x o =1 . 

So when c, s, and V are chosen in this way, we call system (|5.2[) holonomy-rescaled 
locally M. N -invariant Ricci flow and denote a solution by g{r) = (u(r), A{t), G(t)). 

Lemma 8. Let g be a metric of the form h5.1\) such that u = 1, A vanishes, and 
G is a harmonic- Einstein metric. Then g is a stationary solution of L5.2\) . 

Proof. The hypotheses on u and A imply that V = at r = 0. The fact that G is 
harmonic-Einstein implies that |c\)G| 2 is constant in space. Hence i |9oG| 2 = su. 
The result follows. □ 



Notice that the sol-geometry manifolds in Example [3] satisfy the hypotheses of 
the lemma. There s is the topological constant s = 4(logA) 2 . Clearly, s > 
whenever Ma has nontrivial holonomy. 



5.2. Linearization at a stationary solution of holonomy-rescaled flow. Let 

(u + v, B , G + h) be a perturbation of a stationary solution g = (1, 0, G) of the type 
considered in Lemma [8] 

Lemma 9. The linearization of holonomy-rescaled locally M. N -invariant Ricci flow 
about a stationary solution (1,0, G) with G harmonic- Einstein is the au- 
tonomous, strictly parabolic system 

d 

— {v B h)=L(v B h) = (L Q v + Fh L X S L 2 /i) , 



where 






(5.7a) 




= COqV — sv, 


(5.7b) 


Fh 


= (d Q G,doh}-(d%G,h), 


(5.7c) 


(UBY 


= dlB i - S -B i , 


(5.7d) 


(L 2 /i)ij 


— ^ohij — G M (dohikdoGij + doGikdohej) + doGikh M doG 



Proof. One may suppose that Jjr| e _ u = v, ^| e _ A = B. and ^| e _ G — h. We 
first compute that 



d_ 

de 



V = -d v dx° + G ll d B j dx l 
9 
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and 



7k 
d_ 

~de 
d_ 

de 



(£vg)oo = COqV, 



e=0 



£ = 



£=0 



(Cvgho = G^B\ 
(£vg) lj = ^(d v){d () G lj ). 



Note that we used (|5.4p to get the last equality. Using the fact that G is harmonic- 
Einstein, we obtain 



d_ 

de 



\d G\ 2 = 2(d G,d h)-2(d%G,h) 



£ = 



where the inner products are taken with respect to G. Observing that 



d_ 

de 



^oGij — dghij — (doGij) 



£=0 

we conclude that 



d_ 

de 



r° ) =d*h ij -\(d v)(d G ij ), 



d_ 

7Ts 



{u-^lGij - (d G lk )G ke (d G je )}} 



£ = 



d h tJ - G ' (d h ik d G ej + d G ikdo h ej ) + d»G lk h" 8 oGej - -(d v)(d oGij ). 



The result follows. 



□ 



5.3. Linear stability of holonomy-rescaled flow. To investigate linear stability 
of (|5.7p . it is convenient to 'cut' the bundle A4a and consider tensor fields depending 
on x° s [0,1]. By our assumption that Ma admits a flat connection, we may 
henceforth assume that the R^-valued 1-form A has trivial holonomy, i.e. that 
A l \ x o =1 = A % \ x o— , while G satisfies (|5 .€>[) for a nontrivial action of A 6 Gl(iV,M). 
Notice that these assumptions are satisfied by sol 3 -twisted torus bundles, of which 
Example [3] is a special case. 

Corresponding to these assumptions, we require v and B to satisfy the boundary 
conditions 



(5.8) 
(5.9) 



U | x o = ;i — w|a,o = o, 



B' 



E° = l 



B l 



respectively. We require h to satisfy the linear compatibility condition 



(5.10) 



lij |a;° = l 



A^Ajhki | a;o = o 



derived from (|5.6p . 

For now, we apply the differential operators in (|5.7[) only to smooth data satis- 
fying these boundary conditions. Below, we will specify larger domains in which 
smooth representatives are dense. With this intention, we now investigate the lin- 
ear stability of (|5.7|) . We will eventually have to consider the entire action of the 
lower-triangular operator L, but we begin with observations about its component 
operators. The first two are self evident. 



STABILITY OF LOCALLY R -INVARIANT SOLUTIONS OF RICCI FLOW 



23 



Lemma 10. Lo is self adjoint with respect to the inner product 

(<P,i>) = / 
Jo 

and is strictly linearly stable if s > 0. 

Lemma 11. Li is self adjoint with respect to the inner product 

(</j,V>) = / <p l ^ j 5 i3 dt 
Jo 

and is strictly linearly stable if s > 0. 

Lemma 12. L2 is self adjoint with respect to the inner product 

Jo Jo 
and is weakly linearly stable. Its null eigenspace is 

{GM : M G gl(-/V,M) commutes with all (?(•)}■ 

Proof. If (p and ip satisfy (|5. 10(1 and are smooth, then doipij \ x a =1 = A^Aj^oVfe^ U°=o 
and doipij\ x "=i = ^iMdoipke\x°=o- So if G satisfies (|5.6j) . then all boundary terms 
arising from integrations by parts vanish. Using this fact, it is easy to verify that 

(L 2 ip,ip) = (y,L 2 ^) • 

Now we write (L2/1, h) = I± — 2I 2 + ^3, where (using primes to denote differen- 
tiation with respect to x° — £ inside an integral) the distinct terms are 



h = J h'^hi, <;■'■(;■' dc 

h= f ti vp G™G qj h M G lk & 1 d£, 
Jo 

[ G ip G pq h qj G' kr G rs h se G ie G :ik d£_. 
Jo 



h = 

An integration by parts shows that 

h=~ WW 2 + 2I 2 . 

The integral -Z3, which may be rewritten as 

h = (d Goh,hod G), 

requires more work. Integrating it by parts produces seven terms: 

h = 2/3 + 2 \\d G o h\\ 2 - 2 (d G o h, d h) - (dgG oh,h), 

Because G satisfies the harmonic-Einstein equation (|5.5p . the final term above may 
be rewritten as 

(dlGoh,h) = \\d Goh\\ 2 . 
Collecting terms and applying Cauchy-Schwarz, we obtain 

(L 2 h, h) = - \\d h\\ 2 - \\d G o hf + 2 (8 G o h, d h) 

< - \\d h\\ 2 - \\d G o h\\ 2 + 2 \\d G o h\\ \\d h\\ < 0. 

One has equality if and only if = G M hij is independent of x°, hence if and only 
if hij = G lk Mf for some M S gl(AT,R) that commutes with all G(x°). □ 
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Remark 5 . An element GM of the nullspacc of L2 will itself satisfy the harmonic- 
Einstein equation (15. 5|) if it is a metric, i.e. if M is invertible. 



Remark 6. For the sol-Gowdy metrics considered in Example [3l it is easy to verify 
that the nullspace of L2 is two-dimensional. 

Now we extend the L 2 inner products above to vectors (y B h) in the natural 
way, i.e. 

((« B h),(a (3 y))=(v,a) + (B,0) + (h, 1 ). 
The lemma below is the main observation of this subsection. Observe that it does 
not directly imply convergence of the nonlinear system (|5.2I) , because L is not self 
adjoint. In the next subsection, after introducing suitable function spaces, we shall 
overcome this difficulty by exhibiting L as a sectorial perturbation of the self-adjoint 
operator 

(5.11) L sa :(v B h) i-> (L v LiB L 2 h) . 

Lemma 13. Assume that G has nontrivial holonomy A £ G\(N, K) but admits a flat 
connection. Then there exists C > depending only on G such that the linearization 
L of holonomy-rescaled locally R N -invariant Ricci flow is weakly linearly stable with 
respect to the DeTurck diffeomorphisms induced by V — V(C) defined in i5.3\) . Its 
nullspace is 

{GM : M e gl(A,M) commutes with all G(-)}. 

Proof. First suppose that h belongs to the null eigenspace of L2 . Then because G 
is harmonic-Einstein and h = GM, one has (d^G, h) = (8qG, doh) . Hence Fh = 
and so 

(L(w B h),(v B h))=-C\\d v\\ 2 ~s\\v\\ 2 -\\d B\\ 2 ~^\\B\\ 2 . 

This is strictly negative unless v and B vanish identically. 

For the general case, write h = h° + h • , where h° belongs to the nullspace of 
L2 and h 1 - belongs to its orthogonal complement. The self-adjoint elliptic operator 
L2 has pure point spectrum with eigenvalues of finite multiplicity. So there exists 
A > depending only on G such that (L^ft.^, h 1 -) < —A H/i^H • Hence 

(L(w B h),(v B h)) =-C\\d v\\ 2 -s\\v\\ 2 + (Fh ± ,v) 

-\\doB\\ 2 -^\\B\\ 2 -X\\h^\\ 2 . 

Now for any h, one integrates by parts and uses the harmonic-Einstein identity 
(|5.5|) to see that 

(Fh,v) = - [ (G',h)v'd£. 
Jo 

There exists k depending only on G such that |(G", h)\ < K \h\ pointwise. Thus one 
estimates that 

(L(v B h),(v B ^))<- s || w || 2 _f ||B|| 2 -^||^|| 2 

C\\d v\\ 2 -k 1^1^11 + ^Wh^ 

The quantity in braces is nonnegative as long as C > k 2 /(2A). The result follows. 

□ 
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5.4. Convergence and stability of holonomy-rescaled flow. 

Proof of Theorem [3J We again follow the proof of Theorem [T] mutatis mutandis, 
except for two critical steps. 

Step 1. The complexification L^ a of the operator defined in (|5 . is strictly 
elliptic and self adjoint with bounded spectrum, hence is sectorial by standard 
Schauder estimates. Observe that L c - L^ a = F c , where F c G C(Y 1 1 2 +a , T,° 2 +a ) 
for all a £ (0,1). Namely, F c is a bounded operator from intermediate spaces 
X,j = (Xo,Xi)# and E# = (Eo,Ei),j to Xo and Eo, respectively, with its bound 
depending only on the stationary solution about which we are linearizing. The fact 
that L c is sectorial then follows from classical perturbation results. For instance, 
see [21] Proposition 2.4.1]. 

Step 4. The only change in this step, once again, is how one controls the 
diffcomorphism tp T . Here the argument is much easier, because the base S 1 is flat. 
By (15. 3|) . the vector fields V that generate the diffeomorphisms ip T satisfy 

< C N ||7r s g(r) - 7 (vr c g(r))|| Xi < C N Ce~^ ||g(0) - g^ . 

It is then not hard to see that the diffeomorphisms ip T and (p T must converge. 
Indeed, this is proved in [131 Lemma 3.5]. The remainder of the argument goes 
through without modification. □ 
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